Magnetic Quantum Critical Point in YBa 2 Cu 3 07_,5? 
Resolving the Correlation Length Controversy 
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A mode coupling calculation which previously explained 
the Mott gap collapse induced in cuprates by electron doping 
is applied to the analogous problem of hole doping. A plateau 
in the <f-space susceptibility is found to inhibit the rate of 
susceptibility divergence as T — > 0. This effect could reconcile 
neutron and NMR measurements of the correlation length in 
this regime, and clarify the issue of proximity to a quantum 
critical point. 

There is currently a controversy as to the magnetic 
correlation length £ in hole doped cuprates, with neu- 
tron scattering [1] finding a T-indepcndcnt value while 
NMR studies [2,3] suggest that £ increases as T decreases. 
While the numerical values are not greatly dissimilar, the 
temperature dependence is important, since a weakly- 
diverging £ - even if ultimately cut off by superconductiv- 
ity - could provide evidence for proximity to a quantum 
phase transition (QPT), as has been suggested to occur 
slightly above optimal doping [4]. Here, a resolution of 
this problem is suggested, which is consistent with the 
theoretical model of the QPT. 

Recent evidence for a QPT in electron-doped cuprates 
[5,6] has generated a lot of theoretical interest [7-10]. 
It is found that at a critical electron doping, there is a 
QPT where (a) the Mott gap collapses, (b) the T=0 Neel 
transition terminates, and (c) the Fermi surface crosses 
over from small pockets to a large barrel. The simplest 
mean field theory [7] can reproduce the transition, if a 
Kanamori-like [11] doping dependence of the Hubbard 
U is assumed. The results have been confirmed by a 
number of calculations [8-10] which incorporate interac- 
tions with spin waves, the main difference being that the 
mean-field gap and Neel transition become a pseudogap 
and crossover temperature T*. [Secondary interactions 
can generate a finite Neel temperature Tn << T* .] 

Extension of the model to the hole-doping case is 
problematic. The model predicts a nearly electron- 
hole symmetric QPT, with pseudogap collapse. In fact, 
the observed pseudogap does follow the predicted dop- 
ing dependence [8], terminating in a QPT [4] which 
is approximately electron-hole symmetric. This sym- 
metry can be seen in the phase diagrams [12] of 
electron-doped Nd2- a ;Ce :E CuC)4± 1 5 (NCCO) and hole- 
doped La2- a; Sr a ;Cu04 (LSCO). Signatures of strong local 
magnetic couplings disappear near the dopings of opti- 
mal superconducting T c ; the corresponding dopings are 
comparable, and the T c values themselves do not greatly 
differ. Near optimal doping the normal-state resistiv- 
ity is linear in temperature [13], suggestive of a QPT. 



Moreover, it was recently suggested [14] that this QPT 
involves a crossover to the large Fermi surface. 

However, the third element of the transition, the T = 
Neel transition, appears to be absent: the anticipated 
correlation length £ divergence as T — > is cut off at a 
surprisingly low doping in neutron measurements [1] . On 
the other hand, the T-dependent £ deduced from NMR 
could suggest that a weak T — > Neel transition is still 
present. The questions then are two: can the NMR and 
neutron results be reconciled, and why is the Neel transi- 
tion so much weaker in hole-doped cuprates? Both ques- 
tions can be answered by a self-consistent renormaliza- 
tion (SCR) theory [8,15]. 

Within the SCR approach, the zero-temperature anti- 
ferromagnetic (AFM) transition is controlled by a renor- 
malized Stoner criterion [8], Uxoq — with U the 
Hubbard U, r\ > 1 a quantum correction, and xoq = 
XoiQ^ — 0) the bare magnetic susceptibility at Q — 
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where / is the Fermi function, and E the electronic dis- 
persion, taken to be of tight-binding form 
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with Cj = cos kid and a is the lattice constant. At fi- 
nite temperatures, there is a Mott-like pseudogap in the 
density of states, which first appears at a temperature 
close to the RPA Neel temperature (for U* = U/rj) - in 
particular, the pseudogap collapses to T = at the zero 
temperature Neel transition. While the doping at which 
this QPT occurs depends sensitively on |t'|, the QPT falls 
at approximately the same doping for both electron and 
hole doping. 

It is here suggested that the striking differences be- 
tween electron and hole doped cuprates can be under- 
stood from the Stoner criterion, specifically, the doping 
dependence of Xoq for q near Q, Fig. 1. For hole doping, 
the susceptibility has the form of a (diamond-shaped) 
plateau centered on Q; for electron-doping the plateau 
width has shrunk nearly to zero. Briefly, as T decreases, 
the Stoner enhancement arises only very near q = Q for 
electron- doped systems; as the peak height diverges the 
width (inverse correlation length) shrinks to zero, and the 
area under the curve remains small, leading to a strong 
transition, £ diverging exponentially in l/T. For hole- 
doping, the Stoner enhancement is spread over the full 
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plateau. This has two effects: first, the susceptibility 
width is pinned at the plateau width, and is insensitive 
to the £ divergence; second, the large area over which the 
susceptibility is growing can lead to a sum-rule satura- 
tion and a very sluggish divergence, £ ~ 1/vT. The 
detailed calculations below confirm this scenario, and 
show that the £ measured from NMR has the expected 
T-dependence. Thus, the SCR model both explains the 
seemingly contradictory neutron and NMR data, and si- 
multaneously explains the electron-hole asymmetry. 
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FIG. 1. (a) Susceptibility xoq for several hole dopings in 
the hot spot regime, at the mean-field Neel tempreature Tjy. 
From lowest to highest curve, dopings x [and Tn] are [4700], 
0.05 [2200], 0.078 [1520], 0.109 [961], 0.127 [701], 0.147 [475], 
0.160 [369], 0.175 [278], 0.188 [222], 0.214 [133], 0.225 [105], 
0.236 [70], 0.241 [55], and 0.246 [40K]. (b) Susceptibility X 0q 
for several hole dopings x > xh near Q, at T = 10K. From 
bottom to top (at Q): x = 0.300, 0.286, 0.274, 0.261, 0.253, 
0.250, 0.247, and 0.246. (c) Susceptibility xoq for several elec- 
tron dopings near the C-point: for dashed lines from bot- 
tom to top: x [T N ] = -0.147 [1661], -0.171 [1304], and -0.196 
[853K]. Solid lines, from top to bottom: x — -0.213 [C-point, 
T N = 199K], and, beyond C-point, T = WOK, x = -0.223, 
-0.246, -0.269, and -0.293. 

The susceptibility plateau is controlled by hot spot 
physics, where a hot spot is defined as a point on the 
Fermi surface (FS) which is separated from another FS 
point by exactly the antiferromagnetic vector Q. For 
any t' ^ 0, the hot spots exist in a finite range of doping 
about half filling; the critical dopings which terminate the 
hot spot regime are called xc (electron doping) and xh 
(hole doping) - the latter coinciding with the van Hove 
singularity (VHS). Between these two dopings, the sus- 
ceptibility is large [16-18], but outside this doping range 



the susceptibility falls rapidly, so for a wide range of cou- 
plings, U, the Stoner criterion will hold when hot spots 
are present, but will fail when the hot spots disappear. 
The approximate electron-hole symmetry of the hot-spot 
regime leads to a corresponding symmetry of the QPTs. 

Figure 1 shows the bare susceptibility as a function 
of q both in and out of the hot spot regime. At each 
doping in the hot spot regime, there is a plateau in q- 
space centered at Q, where Re(xo q ) is nearly constant, 
Fig. la. This plateau is highly asymmetric between elec- 
tron and hole doping: the width of this plateau q c shrinks 
to zero at xc, Fig- 1c, but actually has its largest value 
at xh- The g-plateau half- width q c is the point where 
the Q + (f c -shifted-FS no longer overlaps the original FS. 
For displacements along the [110] direction, q c is found 
from 

sin — = \/— r , (3 

with t = 2t'/t, while along the [100] direction, q c — 
2sin~ 1 (—fi/2t). The plateau is a Fermi surface caliper: 
Eq. 3 represents the distance in g-space between the 
nodal point and the zone diagonal. The cutoffs are plot- 
ted in Fig. 2. The experimental data [19,1] will be dis- 
cussed below. 

For doping off of the plateaus, either for x > xh, 
Fig. lb, or x < x c , Fig. lc, there are two independent 
QPTs. First, the magnetization vector changes from 
commensurate (q = Q) or weakly incommensurate on the 
plateau to strongly incommensurate off of the plateau, 
due to Kohn anomalies [20,21]. However, while this in- 
volves a transition between two different forms of mag- 
netic order, there is a second transition. Due to the rapid 
falloff of the magnitude of xo> there is a Slater- type tran- 
sition to a nonmagnetic state. It is this latter transition 
which matches experimental observations for the cuprate 
QPTs. The falloff is so sudden that the strongly incom- 
mensurate susceptibilities have not yet been seen. 
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FIG. 2. Widths of susceptibility plateaus near Q, as a func- 
tion of chemical potential {j,. Data from Ref [19] (circles). 

Figure 2 compares the theoretical plateau widths near 
the -ff-point with the measured [19,1] plateau widths 
in YBa2Cu307_5 in the normal state, along the [In- 
direction. [Note that the broad feature analyzed in 
Ref. [19] is distinct from the resonance peak.] For low 



2 



doping ({j, > — 0.28ey) the widths are only weakly in- 
dependent, and are measured in the low-w limit. For 
higher doping there is a spin gap at the lowest frequen- 
cies, but since this may be a sign of a competing insta- 
bility, the magnetic linewidth can be estimated at the 
lowest frequencies above the gap. Particularly at high 
doping, the agreement with the calculated plateau width 
is quite good: i.e., neutron scattering is measuring the 
plateau width, which is not indicative of the correlation 
length. The 'deviations' near half filling (/i = — 0.2eP) 
are in fact the expected SCR behavior, with the system 
developing long-range Neel order, £ ~ l/(Ag) — > oo. In 
turn, the large values of q c can explain why T/v — > at 
such a low hole-doping: the dotted line in Fig. 2 shows 
the criterion £/a = 100, which is the threshold for a finite 
Neel temperature in the electron-doped cuprates [22]. 
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FIG. 3. (a) Normalized susceptibility as the QPT is ap- 
proached, for x = 0.1, T = WOK . From broadest to narrow- 
est, the susceptibility denominator is So = 1, 0.51, 0.27, 0.145, 
0.023, 0.011, 0.0048, 0.0011, and 0.00027. (b) Peak width Aq 
vs normalized correlation length £ = 1/vfc. 

The relation between linewidth Ag and correlation 
length £ can be readily demonstrated. In a self-consistent 
renormalization (SCR) calculation [8,15], the RPA sus- 
ceptibility is renormalized to xq = Xoq/^j with S = 
1 — U Xoq + A, an d A is a fluctuation-induced correction 
which keeps 6 > for T > 0. For q = Q + (f close to Q, 

SocC 2 + q' 2 - (4) 

Thus if Eq, 4 holds over a wide g-range, the peak width 



is a measure of the magnetic correlation length. In the 
presence of the susceptibility plateau, this is no longer 
the case, Fig. 3. Figure 3a shows the SCR susceptibil- 
ity (normalized to its peak value) at a typical hole dop- 
ing, for several values of Sq = S(q' = 0). In Fig. 3b, 
the resulting half-width Ag is plotted as a function of a 
normalized £ = 1/\/Sq. It can be seen that for an ex- 
tended range of £ the linewidth is unrelated to £ and 
measures the plateau width, as expected from Fig. 2. In 
the present calculation, the magnetization becomes in- 
commensurate, causing Ag to start decreasing again and 
— > as £ — > oo. However, this is due to a giant amplifica- 
tion of fine structure on the plateau close to the ordering 
transition; in reality such fine structure may be washed 
out by fluctuations. 

The presence of the susceptibility plateau can affect 
the T-dependence of £. In two-dimensions, fluctuations 
prevent a finite temperature Neel transition, but in the 
renormalized classical regime the susceptibility is ex- 
pected to diverge exponentially as T — > 0, leading to a 
similar divergence of £. However, the susceptibility must 
satisfy a sum rule, and if the peak width stops decreas- 
ing (due to the plateau), this sum rule can be saturated, 
leading to a weaker temperature dependence for £. The 
fluctuation-dissipation theorem can be written as [15,23] 



< M >- 



duo f d 2 i 

ix v > J (2tt 



d 2 q .Cx 



—)Imx{q,u) (5) 



where < M 2 > is the mean square local amplitude of 
nearest neighbor spin fluctuations and n is the Bose func- 
tion. The integral can be separated into a part over 
the plateau and a background term from regions far 
from [24] Q — (jr,Tr), and approximately evaluated as 
A/T — B = q 2 /So where q c is the plateau width, or 



(6) 



which holds at low temperatures, T << a/b. Figure 4 
plots Eq. 6 for the NMR derived coherence lengths [3] for 
optimally doped T c = 90K, and underdoped T c = 66K 
YBa 2 Cu 3 07_,5 (YBCO). It can be seen that both curves 
follow Eq. 6 at low T. Hence, the plateaus explain the 
weak divergence of the correlation length in hole doped 
cuprates. (Note that if there were no plateaus, q c — * 0, 
there would be no constraint on £ - the width would de- 
crease as the peak height increased - and the exponential 
divergence would be recovered, as in the electron-doped 
cuprates.) Finally, since a is proportional to the square of 
the staggered magnetization, the decrease of the a coeffi- 
cient with doping may signal the proximity of the system 
to a magnetic QPT. 
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FIG. 4. Measured [3] magnetic coherence lengths in 
YBCO, plotted according to Eq. 6. Full circles for under- 
doped, open circles for optimally doped samples. Straight 
lines are fits to Eq. 6 with parameters a, b = 8400K, 64.7 
(960K, 3.93) for T c = 66K (90K). 

Thus the present model can explain the high-T prop- 
erties of the pseudogap in the hole-doped cuprates. At 
lower temperatures the hole-doped cuprates are also sus- 
ceptible to a number of competing instabilities, most 
notably superconductivity, but also including stripe 
physics. Certainly, the broad g-plateau is conducive to 
incommensurate modulations, Fig. 3. While these com- 
peting orders will complicate detailed calculations, they 
are unlikely to significantly modify the picture of the 
Mott gap collapse presented here. 

In conclusion, the broad q-plateaus near the H-point 
can reconcile the neutron scattering and NMR measure- 
ments of the correlation length - NMR should see a di- 
vergence, albeit weak £ ~ T^ 1 / 2 , as T — > 0. The neutron 
linewidth measurement should not see such a divergence, 
but the divergence should be reflected in neutron mea- 
surements as a I /T growth of the peak intensity. A reso- 
lution of the correlation length problem along these lines 
would mean that the same Mott gap collapse can explain 
the QPTs for both electron and hole doped cuprates - 
indeed the disappearence of hot spots provides a natu- 
ral phase boundary for QPTs. Superconductivity near an 
AFM or ferromagnetic QPT has recently been observed 
in a number of systems [25] . 
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